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Abstract. The "twisted Mellin transform" is a slightly modified version of the usual classical 
Mellin transform on L-^([0, oo)). In this short note we investigate some of its basic properties. 
From the point of view of combinatorics one of its most interesting properties is that it intertwines 
the differential operator, df/dx, with its finite difference analogue, V/ = f{x)—f{x—l). From the 
point of view of analysis one of its most important properties is that it describes the asymptotics 
of one dimensional quantum states in Bargmann quantization. 



1. Introduction 
The standard Mellin Transform is defined by the formula 



(1.1) Mf{s) = / fix)x'-'dx. 

Jo 

This paper has as its topic the following "twisted" version of ()1.1|) : 

/o°° f{x)x^e~^ dx 



(1.2) 



Mf{s) 



Jo 



dx 



This transform has a number of remarkable properties, the most remarkable perhaps being that it 
intertwines the standard differential operator, and the finite difference analogue of • 



dx ' 



(1.3) 



^f{x) = f{x)-f{x-l). 



By a theorem of MuUin and Rota it is known that there exists an invertible operator intertwining 
the "umbral" calculi generated by ^ and V; but, as far as we know the explicit expression (|1.2p 
for this intertwincr is new. 

The main topic of this paper is an asymptotic formula for (|1.2p when the function, /, in the 
integrand is a symbol of degree fc, i.e. has the property that for any r e N there exists a constant 
Cr such that 



(1.4) 



dx'' 



(x) 



< Cr 



More explicitly, we will show that for such functions 

(1.5) A^/(x)~^/«(x)g.(x), 
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where /*•''•' (x) = "E^fi^)^ ^^'^ 9r{x) is a polynomial of degree [r/2] given by a simple recursion 
formula. In some joint work with V. Guillemin, now in progress, we will use this formula to obtain 
results about the spectral density functions of toric varieties. 

A few words about the organization of this paper. In section 2 we will prove some elementary 
facts about the domain and range of M , derive a twisted version of the standard inversion formula 
for the Mellin transform, prove that M has the intertwining property that we described above and 
compile a table of twisted Mellin transforms for most of the standard elementary functions. In 
section 3 we will use steepest descent techniques to derive (jl.Sp and give two rather different recipes 
for computing the 5r's, one analytic and one combinatorial. (By comparing these two recipes we 
obtain some curious combinatorial identities for the Stirling numbers of the first kind.) 

We would like to thank Richard Stanley for a number of helpful comments on the umbral cal- 
culus, Stirling numbers, and his suggestion on combinatorial properties of the sequence of functions 

2. The Twisted Mellin transform 

Let C be equipped with the Bargmann measure /i = e^'^' dzdz. Given a function / e C°°(C), 
one would like to study the asymptotics of the spectral measure 

nif) = TrK-M^^fc), 

associated with the quantum eigenstate, z^ ^ as fc — > oo, where tt/c is the orthogonal projection from 
i2(C, ^) onto the one dimensional subspace spanned by z^ , and Mf is the operator "multiplication 
by /". By averaging with respect to the T^-action, we can assume / € C°°(C)^\ i.e. 

where r = \z\ is the modulus of complex number z. 
For fc e N, one has 

{fz\z^), fir^y^+^e-^' dr 



Tkif) 



/o°° f{x)x'^e~'-^ dx 



x^e~^ dx 

So the asymptotic properties of Tk{f) can be deduced from asymptotic properties of the twisted 
Mellin transform 

(2.1) MHs) = f*i^) = S^fi^p^ 

Jp x-^e ^ dx 

where / G C°°(R+) and s > 0. Note that the denominator is the Gamma function T{s + 1), while 
the numerator is just the Mellin transform of the function xe~^ f{x). 

For this integral to converge, we will assume that / is of polynomial growth, i.e. 

(2.2) |/(a;)| < Cx^ for some N. 
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Some basic properties of the transform are the foUowing 

Proposition 2.1. Suppose a, 5 G M, c > 0, n G N, / zs a function of polynomial growth, then 

(1) Forgix)^xyix), 

(2.3) Mg{s) = lil±^±}lMf{s + a), 

and for g(x) = e^'^'^ f(x), 

(2.4) Mg{s)^{c+l)-'-^Mfc{s), 
where fc{x) is the dilation, fc{x) = /(^^)- 

(2) Forg{x) = £{x), 

(2.5) Mg{s) = WMfis) := Mf{s) - Mf{s - 1), 
and more generally, for any n € N and g{x) — /'"^(x), 

(2.6) Mg{s) = V-{Mf){s) = f^(-l)^ (")-^/(^ " 

(3) For g{x) = f{x) In a;, 

(2.7) ^Mfis)^Mgis)-Mfisf-^. 

(4) forg(x)=/;/(i) di, 

W-i 

(2.8) M<7(s)= ^ X/(s-*) + Xg(s-[s]). 

1=0 

/n particular, 

n 

(2.9) M<?(n) = ^X/(z); 

i=0 



Proof. The assertion (1) is obvious. 

To prove (2), we note that for g{x) = f'{x), 

Jp x^e"^ ax 

fo° f{^){^'' ~ sx^~^)e~'' dx 
x^e~~^ dx 

f{x)x'^e~^ dx f{x)x'^~^e~^ dx 

x^c^^ dx x^~^e~^ dx 

= Mfis)-Mf{s-l). 
The property (|2.6p is easily deduced from (|2.5p by induction, and (3) is a direct computation. 
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To prove (4), we note that by integration by parts, 

(2.10) Mg{s)=Mf{s)+Mg{s-l), 

which imphes As for (P^ . this foUows from the obvious fact Mg{0) = Mf{0). □ 

From the definition its easy to see that the twisted Melhn transform is smooth, i.e. it transform 
a smooth function to a smooth function. Moreover, it transforms a function which is of polynomial 
growth of degree iV to a function which is of polynomial growth of degree N, and Schwartz functions 
to Schwartz functions: 

Proposition 2.2. (1) Suppose \f{x)\ < Cx^ , then \Mf{s)\ < C's^. 
(2) A4 maps Schwartz functions to Schwartz functions. 

Proof. (1) This comes from the definition: 

rCx^^-e- dx _ rjs + N + l) 
'■^^^'^1- T{s + 1) T{s + 1) • 

(2) Suppose / is a Schwartz function, i.e. for any a, (3, there is a constant Ca,p such that 
supjx"9'^/(x)| <C„,0. 

For /? = 0, \x°'f{x)\ < C implies |s"A^/(s)| < C . 

For /3 = 1, we apply ([27]) and the above result to get \s°'^^Mf{s)\ < Cq. 

For /3 > 1, let ip{s) = r'(s)/r(s). Then by repeated applications of (|2.7p one can see that 
-^Mf{s) is a linear combination of the functions Mgi{s)')p'''^'^ {s + 1), where gi{x) — f{x){hixy 
and 

(2.11) ^M(, + 1)^_V,(. + 1) 



is the polygamma function, which is bounded for each m, as is clear from its integral representation: 
Thus by induction we easily deduce that \s°'d'^Aif{s)\ < Ca,p. □ 



4-m -(s+l)t 



</ — -dt = C(™ + l)r(m+ 1). 

1 — e ^ 



Remark. Since the twisted Mellin transform transforms a Schwartz function to a Schwartz function, 
we can define the twisted Mellin transform on tempered distributions by duality. 

We will next compute the twisted Mellin transform for some elementary functions such as 
polynomials, exponentials and trigonometric functions. 

(a) For f{x) = x" , 

(2.12) X/(s) = r(s + a+l)/r(s + l). 
In particular, if f{x) = x", n a positive integer, then 

(2.13) Mf{s) = si"! .= (s + + 2) . . . (s + 7i). 
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Thus the twisted Melhn transform of a polynomial of degree n is again a polynomial of degree n. 

(b) Suppose a > 1, then for f{x) — , 

(2.14) X/(s) = (lna + l)-l-^ 
More generally, if f{x) — x^a^^, then 

(2.15) Mfis) = (Ina + l)-i-''-T(s + b + 1)/T{s + 1). 

(c) For fix) = j^, 

(2.16) Mf{s)^C{s + l), 
and as a corollary, for the Todd function f{x) = j^p^, 

(2.17) M/(,s) = (s + l)C(s + 2). 

(d) For f{x) — Ina:, one gets from p.7|) 

and in general, for f{x) — (Inx)", 

(e) For the trigonometric functions f{x) = smx and g{x) = cos a;, 

Mf(s) = — ^ sm^ 

^ ^ (\/2)«+i 4 

(2.20) 

/ ^ 1 (s + l)7r 
7\4g(s) = — ^ cos ^ 

' (%/2)«+i 4 

[Prooi. Let /i(a;) = e", then Mh{s) = ■^Y^jprr, which gives ^TMl-) 
Similarly for ,f{x) — sin(aa:) and g{x) ~ cos(ax), 

Mf{s) = (1 + a^)^'* sin(s arctana), 

(2.21) 

Mg{s) = {1 + a'^) * cos(s arctana). 

Some concluding remarks: 

(1) From the inversion formula for the Mellin transform, we obtain an inversion formula for 
the twisted McUin transform: 

(2.22) f{x) = e'' r{s + l)Mf{s)x-'-Us. 

J c—ooi 

and from the Parseval formula (c.f. ^) a "Parsevel-like" formula for Ai 



(2.23) / fix)gix)x'e-'^ dx = - Mf{l - s)Mg{s)^. ds. 
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(2) Letting (x)^") = x{x + 1) ■ ■ ■ {x + n - 1), (j^l^ becomes Mf{s - 1) = (s)("\ where 
f{x) — x". Also expanding s^"! in terms of s", we get 



(2.24) s["l = (s + l)(.s + 2)---(s + n) ==^c(n + l,/c + l)s'= 

k=0 



where c(n, fc) is the signless Stirling number of first kind, c.f. [?]. Note that both {a;"} and {x'"^} 
are a basis of the polynomial ring, thus is a bijection from the polynomial ring to itself. 

(3) Formula (|2.5p tells us that A4 conjugates the differential operator, to the backward 
difference operator p.3p . fn combinatorics both ^ and the backward difference operator are 
"delta" operators, with the functions x" and x*^"^ as their sequence of basic polynomials. Thus by 
a theorem of R.Mullin and G-C.Rota ([2]), the map T : /(x) A4f{s — 1) is invertible and the 
map 5* I— > TST^^ an automorphism of the algebra of shift-invariant operators onto the algebra 
of polynomials. Moreover, T maps every sequence of basic polynomials into a sequence of basic 
polynomials. Such an operator is called an umhral operator in the umbral calculus. 

(4) If we replace the Bargmann measure, /i, by the generalized Bargmann measure /Xq, = 
g-"!^! dzdz, then we are naturally led, by the argument at the beginning of this section, to 
studying the "a-twisted Mellin transform" 

(2 25) M f(s) r/(^)^^^'"^ dx 

All the properties in Proposition 2.1 can be easily generalized to Aia- Moreover, it is easy to see 
that 

(2.26) Mo^fcis) = M^/Jis), 

where /c(x) — f{cx), so the a-twisted Mellin transform of a function can easily be computed from 
the twisted Mellin transform. 

(5) The twisted Mellin transform can also be defined in higher dimension in a similar way. 
For the function /(ri, • • • , r^) = fi{ri) ■ ■ ■ fdi^d), its twisted Mellin transform is just the product 
of the twisted Mellin transform of /i, • • • , fd- 



3. The Asymptotic Expansion. 
We can rewrite the twisted Mellin transform as 

r.n Uf(^ /o°°/(x)e-^°g--dx 

(3-1) Mm = ^^oo^,,,^,_,^^ ■ 

For the phase function <y9(x, s) = slogx — x, we have 
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thus the function (ps{x) — ip{x,s) has a unique critical point at x = s. Moreover, this is a global 
maximum of (p{x, s), since 

lim 1^9(3:, s) = — cxD, 



and 



X — ' + 00 



5V 



Hence if / is a symbol, we can apply the method of steepest descent to both denominator and 
numerator to get 

(3.2) Mf{s)^Y.9kis)f^'\^)- 



To compute the functions gk{s) consider the Taylor expansion of /, 

00 

(3.3) 

r=0 

Applying M to (3.3) with s fixed we get 
(3.4) 
where 



r=0 



(3.5) 
For r < 5 small 



00 ^ 



x'^e ^ dx 



his 
his 
his 



= /l(s) = l, 

= 2 + s, 

= 6 + 5s, 

= 24+ 26s + 3s^ 



= 120+ 154s + 35s^ 

which suggests that /.r(s) is a polynomial of degree [r/2] in s. We will give two separate proofs of 
this fact, one combinatorial and one analytic. 
Proof 1. Putting (PTM)) into we get 



i=0 k=0 



j=0 \i=j 



cii + l,i- j + 1) s'' ^. 



On the other hand, by the definition of the Stirling number, 
(3.6) c(n + 1, 71 + 1 — fc) = c(n, n — k) + nc{n, 71 — (fc — 1)). 
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and from this recurrence relation we will show: 

Lemma 3.1. There are constants Cij, depending only on such that 
(3.7) c{i + l,i+l-3)=Y,Ci,j {i)u 

where (i)i ^ i{i — 1) ■ ■ ■ {i — I + 1) is the falling factorial. 

Proof. This is true for j = 0, since c{n + 1, rt + 1) = 1. Notice that 

P'-"(C)H'r)--(:))-<::;)^i^<-'- 

Now use induction and the recurrence relation (13.61). □ 



Now suppose 2j < r, then the coefhcients of s^ ^ in fr{s) is 

i=j ^ ^ i=j ^ ^ l=j 

i-l 



l=j i=l 



=0, 

which proves that fr is a polynomial of degree [r/2]. 

Proof 2. First we derive a recurrence relation for fr{s). Using 

d X — s 

— [s log X — X) = 

dx X 

and integration by parts we get 



r(s + l)/,(s) = - / e"i" — .Nr-i a 
Jo 

d 



x(x — sY —(s\nx — x)dx 
dx 











r 



e'''"^-^ — (a;(x-s)'^-i) dx 
dx 

dx 

x"e~^(x - sy-^ dx + (r - l)s / x^er'^ix - sf-^ dx, 

"'0 



(3.8) fris) - r/,_i(.s) + (r - l)s/,_2(s). 
Moreover, we can compute the initial conditions directly 

(3.9) /i(s) = /o(s) = l. 
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Remark. The recurrence relation (|3.8p also follows easily from (|2.3p and (|2.5p . In fact, if we denote 
hr{x) = (x — sY, then fr{s) = Mhr{s), and thus 

rfr-iis) = MKis) - Mhr{s - 1) 

= fr{s) ~ {M{xK-i){s - 1) - sMK-i(s - 1)) 

= fr{s) - s{Mhr-l{s) - Mhr-l{s - 1)) 
= /r(s) - sir - l)A^/lr_2(s). 

From p.Sp . (|3.9p and induction, it follows again that /r(s) is a polynomial of degree [r/2]. 
Thus coming back to (3.4) we have proved 

Theorem 3.2. For any symbolic function f , we have 

(3.10) A4/(s)^5]i^/W(s)/,(5), 

r 

where fr{s) is the polynomial of integer coefficients of degree [r/2] given by {51 

The polynomials /r(s) have many interesting combinatorial properties: 
(1) Since fr{s) is a polynomial of degree [r/2], we can write 

[r/2] 

(3.11) fr{s) = ^ ar,is\ 

i=0 

the coefficients satisfying the recurrence relation 

(3.12) ar,i = rar-i,i + (r - l)a,._2,i-i 
and initial conditions 

0^,0 ==7-!, a2fc,fc = (2A: - 1)!!, 

which implies 



,,1 1 1 

ar,i =f ' - + 



ar,2 =rl 



r r-1 2/ ' 

{r-l)ar-2,i , (r-2)ar-3,i , , 3a2,i 



and in general 
(3.13) ar,k = rl 



r! (r-iy. 4! 

{r - l)ar-2,k-i , (?' - 2)ar-3,fc-i , , (2fc - l)a2fe-2,fe- 



(r-1)! (2fc)! 



(2) The coefficients, a^^i, of /r(s), are exactly those appeared as coefficients of polynomials 
used for exponential generating functions for diagonals of unsigned Stirling numbers of the first 
kind. More precisely, for fixed k, the exponential generating function for the sequence {c{n + 1, n + 
1 — fc)}„>o is given by (c.f. sequence Al 12486 in "The On-Line Encyclopedia of Integer Sequences") 

2k 



n=0 ■ n=fe ^ ' 
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(3) The sequence of functions /^'s have a pretty simple exponential generating function: 

oo ^ oo oc 



^ oo oo 1 / \ 

r! ^ — ' ^ — ' " r\ \i 

r—O i—0 r—i ^ 



\i=0 / \r=i ^ ' 



1+s ■ 



(4) From the generating function above we get a combinatorial interpreting of /r(s) for integers 
s: r!/r(s) is the number of r x r N-matrices with every row and column sum equal to 3 + 2s and 
with at most 2 nonzero entries in every row. (c.f. Exercise 5.62 of 0]). 

(5) There are also other combinatorial interpreting for small value of s. For example, the 
sequence /r(l) count permutations w of {1, 2, • • • , r + 1} such that w{i + 1)7^ w{i) + 1 (c.f. the 
sequence A000255 of "On-line Encyclopedia of Integer Sequences"). For s = 2, we have 



where Dy. is the set of all r x r matrices of ±l's. (c.f. Exercise 5.64(b) of '4^). 

We will conclude by deriving a slight variant of the asymptotic expansion above, which will 
be needed for the application in T]. Given a symbolic function /, consider the integral 



(3.14) 



/Q°°/(x)a;^^e-^" dx 



as N ^ 00. By definition, this is just the "TV-twisted Mellin transform" A4Nf{Ns), which, 
according to fTZE^ . equals MfNiNs), where /at (a;) = f{x/N). Thus by Theorem E21 

(3.15) ^^(/)(s)^^(^l)'/W(s)g..(iVs). 

Note that since gk{x) is a polynomial of degree [k/2], the above formula does give us an asymptotic 
expansion. In particular, we have 

AN{f)is) = f{s) + 1 (/'(s) + /"(s)|) + (^/"(s) + /"'(s)| + f^'\sf-^^ + 0{N~'). 
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